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Zero-dimensional dissipative action possesses non-trivial minima known as Korshunov instantons. 
They have been known so far only for imaginary time representation that is limited to equilibrium 
systems. In this work we reconstruct and generalise Korshunov instantons using real-time Keldysh 
approach. This allows us to formulate the dissipative action theory for generic non-equilibrium 
conditions. Possible applications of the theory to transport in strongly biased quantum dots are 
discussed. 
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I. INTRODUCTION 

Quantum systems coupled to an environment have 
been extensively studied for many years following the 
seminal work by R. P. Feynman and F. L. Vernon. More 
recently a large variety of theoretical models that corre¬ 
spond to different physical realisations of dissipative sys¬ 
tems have been proposed. These models include shunted 
Josephon junction, tunnel junction,' granular array, 1 
Luttinger liquid with a static impurity, ’ and open quan¬ 
tum dot. ’ 

Despite such a diversity the behaviour of these systems 
at low energies turned out to be similar. Indeed, at long 
times the systems become essentially zero dimensional 
and can be described using a collective degree of freedom 
that is often called the phase. In a quantum dot, for ex¬ 
ample, the phase can be viewed as a conjugate variable 
to the electron charge on the dot. The phase dynam¬ 
ics is non-local in time and is governed by Ambegaokar- 
Eckern-Schon type of theory' 1 (AES) that is sometimes 
referred to as dissipative action. 

Linearisation of the dissipative theory reproduces the 
so-called Caldeira-Leggett model : that has a number of 
applications. In particular, the model is instrumental 
for analysing dephasing rate in q-bits and dissipative 
decay rate of metastable states as well as for calculating 
other related quantities. It is worth mentioning that 
similar models were rigorously derived in the context of 
disordered low-dimensional systems as an effective theory 
for low-energy virtual fluctuations. 

It has been realised already in Ref. 14 that periodicity 
of the dissipative action with respect to the phase (which 
originates in the discrete nature of electron charge for 
the quantum dot model) may have some important con¬ 
sequences that are absent in the Caldeira-Leggett model. 
In field theory such non-perturbative effects are often 
accounted for by the so-called instanton solutions. For 
dissipative action those are Korshunov instantons that 
correspond to trajectories of the phase <f> (a zero dimen¬ 
sional bosonic field) labeled by a winding number W. 
This number shows how many times the corresponding 
trajectory circles around the origin. 


Instantons and their contribution to observable quanti¬ 
ties have been thoroughly studied within Matsubara for¬ 
malism that refers to quantum action which is taken on 
the imaginary time contour. In particular, the analy¬ 
sis of the action on the saddle point trajectory has been 
performed in Ref. 15, while the corresponding ’tHooft 
determinant has been found in Refs. 17,18. Later the 
two-loop renormalisation group approach to the prob¬ 
lem has been developed : and the analogy with topo¬ 
logical Pruisken instantons, which appear in topological 
sigma model with broken time reversal symmetry, has 
been established. 1 . The equivalence of the AES model to 
the so-called “paperclip” model with 6 topological term 
is demonstrated in Ref. 21. Recent progress in the field 
is summarised in Ref. 22. 

Instantons distinctly manifest themselves in transport 
and thermodynamic measurements due to their topolog¬ 
ical nature. Instanton solutions in open quantum dots 
describe the emergence of Coulomb blockade s , that man¬ 
ifest itself through weak modulation of the conductance 
with a gate voltage. These are remains of Coulomb block¬ 
ade that dominates transport in closed dots, where con¬ 
ductance vanishes except for a set of narrow Coulomb 
peaks '. In an open dot limit, the Coulomb blockade ef¬ 
fect is strongly suppressed, and results only in oscillations 
of the conductance with a gate voltage. In this case the 
instantons with a certain winding number are responsi¬ 
ble for the corresponding harmonics of conductance os¬ 
cillation. This effect was observed in recent experiments 
with controlled coupling between quantum dot and the 
leads. 

Even though the physics of Coulomb blockade at equi¬ 
librium is fairly well understood its modification under 
generic non-equilibrium conditions remains largely unex¬ 
plored (see, however, Ref. 22). Apart from clear experi¬ 
mental motivation there also exists a substantial theoret¬ 
ical interest to the problem. Non-perturbative analysis 
applied to systems which range from simple quantum me¬ 
chanics to advanced field theory models is equivalent to 
finding “classical trajectories” that minimise correspond¬ 
ing imaginary time action. In contrast, non-equilibrium 
conditions require the use of real time Keldysh formalism 
that may not reconcile with the non-perturbative equi- 


2 


Vg 



FIG. 1: Schematic illustration of the setup: quantum dot, 
source, drain and gate electrodes. 

librium approach. 

Some progress in this direction has been achieved 
within a phenomenological framework utilising quantum 
Langevin equation. While this approach captures many 
qualitative properties of the system it can be justified 
only in the topologically trivial sector while the consid¬ 
eration of instanton contributions requires more careful 
analysis. 

In this paper a non-perturbative approach based on 
Keldysh formalism is developed for a generic zero¬ 
dimensional non-linear theory. A non-equilibrium ex¬ 
tension of the dissipative action for a biased quantum 
dot is constructed and the corresponding integral saddle- 
point equations are derived and solved explicitly in terms 
of hypergeometric functions. For equilibrium conditions 
known results are reproduced. 

In the next Section we shell focus on the formulation 
of dissipative action in real time using an open quantum 
dot as a prime example. 

II. DISSIPATIVE ACTION IN REAL TIME 

Let us focus on a system consisting of a quantum dot 
that is schematically represented in Fig. 1. The electron- 
electron interactions on the dot are described by the 
Hamiltonian H c = E c (n — n g ) 2 , where h is the electron 
density operator and n g is a background density set by 
the potential of the gate. The characteristic charging en¬ 
ergy E c = e 2 / 2C is determined by the capacitance C of 
the dot and by the electron charge e. We assume that 
the dot is metallic, with a single particle level spacing A 
that is smaller than all other energy scales in the problem. 
Each of the two leads coupled to the dot are assumed to 
be kept at thermodynamic equilibrium. The coupling is 
fully characterised by a unitary scattering matrix S(E ), 
which defines all transport characteristics of the dot at 
an energy E in the absence of interactions. 


Coulomb interaction on the dot is decoupled by 
Hubbard-Stratonovich transformation. As a result elec¬ 
trons on the dot are subject to the confining poten¬ 
tials, disorder and time dependent auxiliary bosonic filed. 
The problem at hand can be, therefore, viewed as the 
Ladauer-type transport problem in the presence of ex¬ 
ternal quantum field. Alternatively, the interaction on 
the dot (that is taken in the form of the charging energy) 
can be placed on the leads. This transforms the prob¬ 
lem to a network of interacting quantum wires (scatter¬ 
ing channels) J connected via a non-interacting quantum 
dot described by the scattering matrix S. 

Quite generally all physical observables in the system 
are given by the fermionic path integral 

{O) = J D'FZM' 0{^]e iA[ *\ (1) 

where the real-time fermonic action on the Keldysh con¬ 
tour A[T] = § K dt (z'F<9 t 'k — H) can be decomposed to 
the free fermion and to the interaction parts 

A[V] = Aof'F] + A C [T]. (2) 

Following the standard route we shall perform the 
Hubbard-Stratonovich transformation to decouple inter¬ 
action term by means of the bosonic field <fi and integrate 
out the fermionic degrees of freedom. The resulting the¬ 
ory is given by the functional integral 

{O) = J V(/)0[(j)\e iA W (3) 

where (f> = dt4> is a potential on the dot, which is uniform 
in space by construction. The resulting action is readily 
decomposed into a sum 

A[(j)\ = A c [(j)\ + Adot^]) (4) 

that consists of the charging contribution 

Ac[4>\ = -x4“ dt (/> 2 (f) -n g (j) dt 4>(t), (5) 

C J X J X 

and the coupling term 

A d ot = -i In 2 [</>]. (6) 

The latter is given by the sum of “vacuum loops” and 
can be generally written as 

Z[4>] = det [1 - / + %]/] (7) 

where n a is an advanced component of the fermionic 
polarisation operator, / is a single particle distribution 
function operator, which is diagonal in energy represen¬ 
tation for stationary systems, and 4> c {4> q ) is the classi¬ 
cal (quantum) component of the bosonic field. Namely, 
4> c = {4>+ +4>-)/2 and 4> q = </>+ — 0_, where </>+(-) stands 
for the bosonic held on the upper (lower) branch of the 
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Keldysh contour. The operator R is given by the poten¬ 
tial dependent single-particle scattering matrix 

i2[0] = 5t[^ + ] < S[0_]. (8) 

In general the operator R is non-local in time domain. 

An essential simplification to the general theory of 
Eqs. (5-7) is achieved in the so-called tunnelling limit, 
which corresponds to a weak coupling between the dot 
and the leads. The dependence of the scattering matrix 
on external potential <fi yields in the tunnelling limit 

« = (A./T). »> 

where r(r') and t(t') are the matrices of reflection and 
transmission amplitudes in channel space for a non¬ 
interacting system. (Here we ignored the dependence of 
transmission and reflection probabilities on energy within 
the window given by E c ). 

All transmission probabilities T„, which are the eigen¬ 
values of the matrix tv (or t'(t')i), are small in the tun¬ 
nelling limit T„ C 1. Thus, it is legitimate to expand the 
logarithm in Eq. (6) to the lowest order with respect to 
t and t!. With the help of the identity IndetT = TV In A 
we obtain in the tunnelling limit 

Adot = jf dt 1 dt 2 $L(ti)a(t 1 ,t 2 )$R{t2), (10) 

where g = gL+ 9 r , and g v = '}2 n T% is the dimensionless 
coupling between the dot and the left (77 = L) and the 
right (77 = R) leads, also known as Caildera-Legget “vis¬ 
cosity” and we focus on the limit of a strongly coupled 
dot (g 1); d> a are the vectors in Keldysh space 



where the index ± denotes the upper (lower) part of the 
Keldysh contour and the upper index T stands for the 
vector transposition. The current through the quantum 
dot, can be written as 

j-_ _9 l9r_ f _ y/2) — n(e + V/2)\v(e)/vQ (12) 
gL+gR J 

Here n(e) is a Fermi-Dirac distribution function, v(e) is 
a tunneling density of states, /r are couplings between 
the dot and the left and the right leads. The tunneling 
density of states can be recast in terms of the bosonic 
fields as 7 


v(e) =uq Re 


dre 11 


(1 -n(r))/e^- ( r)-i<M0)^ 



0 )^ 


(13) 


where vq is the density of states for non-interacting elec¬ 
trons. The problem is thus reduced to the calculation of 
the tunneling density of states for this system, both in 
and out of equilibrium. 



FIG. 2: Schematic illustration of the conformal map w in the 
time complex plane, which maps the Keldysh contour K onto 
the contour C. 


A. thermal equilibrium 


In thermal equilibrium the matrix integration kernel in 
Eq. (10) becomes the function of a time difference only, 
a = a(t± — t 2 ), where 


a(t) 


( «++ a + _\ 

V“-+ a —J 


(14) 


is parameterised by the following functions of time 

rj ->2 

ci!++(i) = a _(t) = i Re 1 


a+-(t) = = - 


sinlr(7rTt + iO). 
iT 2 


sinh' (irTt — iO) 


(15a) 

(15b) 


where T is the temperature. 

Analytic properties of the kernel a(t) in Eq. (15) are 
dictated by the shape of the Keldysh contour I \, which 
goes around a cut directed along the real time axis. One 
may reformulate the problem by taking advantage of a 
conformal map of the time complex plane onto the com¬ 
plex plane of the parameter w 21,28 


w = tanh(7rT£). 


(16) 


Under such transformation the strip in the complex plane 
—i/4T < Imt < i/4T is mapped onto the interior of the 
unit circle |ru| < 1 as shown in Fig. 2. Consequently, 
the Keldysh contour K is transformed into the contour 
C that surrounds a cut w £ (—1,1) in the complex w 
plane. 

In terms of the variable w the theory is described by 
the action 

Adot =g / £ J w')$ R (w'). (17) 


The kernel of the integral transform a(w) has the follow¬ 
ing components 


a+.f- = a _= 


1 


+ 


1 


a .|_ (w) = a _(- (—w) = 


(■w + iS ) 2 {w — id) 2 

i 


(w + id) 2 ’ 


(18a) 

(18b) 
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where we introduce a finite shift 6 which implies that 
the corresponding pole is situated outside the countur C. 
The functions z±(w) = z(w±iO) are related to a complex 
function z(w) that is analytic everywere in the complex 
w-plane with the exception of the cut w £ (—1,1). 

One can note that taking the limit of zero temperature 
in the kernel a given by Eq. (15) would lead to the same 
expressions as in Eq. (18) but with the cut extending over 
the entire real axis. Thus the only difference between zero 
and finite temperature analysis is encoded in terms of the 
variable w in analytical properties of the function z(w). 

The action (17) gives rise to the following saddle-point 
equations on the analytic function z(w), 


* dw' 

’z±{w) 

z(w ')' 

z(w') 

z±(w)_ 

1c 


1 1 

(w — w') 2 (w — w' ± iS) 2 




Similarly the instanton solution, which corresponds to 
a winding number W, is given by 


z w (w) 


w 


n 


W-& 

l -£iw’ 


(23) 


with the real parameters £,;, £, £ (—1,1). Substituting 
Eq. (23) in to the action (17) one finds 

A dot [z w ] = igW/2, (24) 

while the charging part of the action is still of the order of 
T/E c corresponding to a generalization of Eq. (22). The 
charging contribution is irrelevant as far as T/gE c <C 
1. By expanding Eq.(12) over the instantons in the low 
voltage limit, one finds a d.c. conductance 


Goc G 


Drude 


l + ^^e g \ w \! 2 cos(2-KWn g ) , 
w 


(25) 


Here the integral over the running variable was closed 
into a countur C , such that the pole at w' = w lays 
inside the contour C, while the pole at w' = w ± iS is 
outside the integration contour. 

Equations (19) were derived by minimizing the action 
Hdot only. Taking into account the charging part of the 
action A c gives rise to the terms which are subleading in 
the parameter T/gE c as shown in the Appendix A). At 
this level of accuracy, i.e. for T/gE c -C 1, the subleading 
terms have to be disregarded. 

Despite the approximations made the saddle-point 
equations (19) are the non-linear integral equations. The 
existing mathematical methods are, however, limited to 
the linear equations only. ’ Fortunately, the restric¬ 
tive analytic structure of Eqs. (19) makes further analysis 
possible. 

It is instructive to take advantage of the analogy with 
the Matsubara representation, which corresponds to the 
contour C taken as the unit circle. Using analytic prop¬ 
erties of the kernel a as described in Appendix B one can 
demonstrate that 


z\ (w) 


w — £ 
1 — £u> 


( 20 ) 


is an instanton trajectory with a winding number equal 
to one, where the parameters £ and £ are independent 
real numbers on the cut £,£ £ (—1,1). The solution of 
Eq. (20) minimizes the action (17). 

The value of the action (17) on the intanton trajectory 
is readily evaluated as 


A dot [z 1 } = ig/2, (21) 


while the charging part of the action evaluated on the 
same trajectory gives 


where 


= e 2 gLgR 
2tt g L + 9r ’ 


(26) 


Note, that the contribution of the instantons with a 
winding number W lead to the appearance of W of the 
d.c. conductance with the gate voltage. Such harmonics 
were recently observed experimentally^ . The fluctua¬ 
tions around the saddle point, account for the zero bias 
anomaly, leading to the renormalziation of the coupling 
strength g —> g + In T/E c , 19,20,31 , and consequently to a 
temperature dependent prefactors in Eq.(25). 

Thus, the construction of real-time Korshunov in- 
stanons using Keldysh formalism at equilibrium is anal¬ 
ogous to that in Matsubara representation once the con¬ 
formal map (16) is employed (see Appendix A for tech¬ 
nical details). The apparent advantage of the Keldysh 
approach is due to the possibility to generalize it for a 
generic non-equilibrium situation. 


B. beyond equilibrium 

To be more specific we shall focus below on the case of a 
quantum dot, which is strongly biased by a source-drain 
voltage V. In this case, the single-particle distribution 
function inside the dot can be parameterized as 

/( e ) = ~/o( e -) + ~/( e +)> (27) 

9 9 

where e_ = e - eVg R /g and e + = e + eVg L lg. 

In terms of the variable w the integral kernel in 
Eq. (17), which depends now separately on two argu¬ 
ments, is given by 


A c [zi] — 2iw + n g- (22) 

The results of Eqs. (21,22) are fully equivalent to those 
obtained in Matsubara representation in the Appendix A. 


onj ( w , w') = aff (w — w') L ^ 

L(x) = 9 ~ l+ / r + 9 -l ^ [ x ™ + x~ iv ] , 


1 — w 1 + w' 

1 + w 1 — w' J 


(28) 

( 29 ) 
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where v = eV/2 -kT and a^\w) is the equilibrium kernel 
set by Eq. (18). The corresponding saddle-point equa¬ 
tions on the analytic function z(w) reads 


i j / r i ^ z ( w ') 

® aw a(w±,w )— 7 —- 

Jc L z±{w) 


a(w', w±) 


Z±{w) 

z(w') 


= 0 , 


(30) 

where z±(w) = z(w ± * 0 ), w± = w±i6, and w £ (— 1 , 1 ). 

We shell look for the solution of Eq.(30) using the fol¬ 
lowing Ansatz 


, (31) 

which is an obvious extension of the equilibrium instan- 
ton of Eq.(20). Motivated by the analytic properties of 
the solution at equilibrium and kernel of integration (15), 
we shell require that the functions F(w) and 1 /F[w) have 
no singularities outside the contour C. That is equiva¬ 
lent to the requirement that all zeroes and singularities 
of F(w) are confined to the interval w £ (—1,1). Under 
such an assumption the saddle point equations can be 
easily solved (see Appendix C for the details). 

A general instanton solution with a winding number 
one is given by an indefinite integral 


z{w) = (1 -££) 


dw 


j (1 — £ w ) 2 

that can be explicitly taken as 


1 — w £ + 1 
1 + w £ — 1 J ’ 


z(w) = K 


+R + 0 9Lgn i-££ 


1 — u>£ 
x ( B v {x) + B-v(x) + 


9 2 1 - £ 

(1 - w){ 1 + £) 
1 — u>£ 


(32) 


(33) 


where K = F( l/£), R is an integration constant and the 
following definitions have been used 


1 — w £ + 1 
X= 1 + £ — 1 ’ 

B v {x) = x lv ( iF 2 ( 1 , iv\ 1 + iv, x) 



(34) 

(35) 


Here, 1 ^ 2 ( 0 , 6 ; c; x) is the hypergeometric function. 

It is worth noticing that the integration constant R can 
be absorbed in the redefinition of the parameters £ and 
K by means of the transformation £' = (£ — i?)/(l — 1?£) 
and K' = K{1 — !?£). By performing this transformation 
we effectively set R = 0 in Eq. (33). The remaining 
parameter K reflects the overall scaling invariance of the 
saddle-point equation (30). Its choice does affect neither 
physical observables nor required analytic properties of 
the function F[w). Hence we can set K = 1 without loss 
of generality. 

Thus, the resulting instanton solution is given by 
Eq. (33) with K = 1 and R = 0. Using known properties 
of the hypergeometric function we can indeed check that 


the solution obtained corresponds to F(w) having ana¬ 
lytic properties required by the construction. For a finite 
value of voltage the instanton solution acquires a branch 
cut in the interval (—1,1). In the limit of equilibrium 
(for gL = 0 , or v = 0 ) the solution reduces to a familiar 
form of Eq. (20). 

Since this interval w £ (—1,1) precisely coincides with 
an image of the Keldysh contour on w plane one may 
calculate the value of the action on the saddle point so¬ 
lution (31), by closing the integrals in the exterior of this 
interval in the complex plane. As a result one gets 

4iot[zi] = j) dlnz(wi) = y . (36) 

Thus, though both, the action and the saddle point so¬ 
lution, are affected by the external bias, the value of the 
action on the saddle point trajectory out of equilibrium 
remains the same. The instanton solutions with a higher 
winding numbers may be constructed in a similar fashion. 


III. CONCLUSIONS AND OUTLOOK 

In this work we have studied the real-time dissipative 
action. We established that the real-time Keldysh and 
imaginary-time Matsubara framework are almost iden¬ 
tical in equilibrium once a conformal map of the time 
plane into wxplane is employed. Topological instantons 
that solve integral saddle-point equations for the Keldysh 
action are constructed. In equilibrium the solutions are 
represented by slightly distorted Mobius maps of w Rie- 
mann sphere onto itself. The index of the map is equal 
to the winding number of instanton. 

The real-time theory is, however, further extended to 
the case of generic non-equilibrium situation. The saddle- 
point equations in non-equilibrium are derived and solved 
analytically to construct the non-equilibrium extension of 
Korshunov instantons. In this paper we have focused in 
particular on the interacting quantum dot that is pushed 
out of equilibrium by applying a large dc bias voltage. We 
demonstrated that, in this case, the saddle point solution 
acquires the branch cut along the entire Keldysh con¬ 
tour. Remarkably, the modified instanton solution of a 
dissipative action survives such a generic non-equilibrium 
condition and the value of the action on the instanton tra¬ 
jectory remains unchanged. It implies, for instance, that 
the conductance oscillations with multiple periods (25) 
emerge when a finite voltage bias is applied. 

Our results may be relevant for a large number of phys¬ 
ical situation that are modelled by the dissipative ac¬ 
tion. To make a more detailed prediction for observable 
quantities such as temperature and voltage dependence 
of electric current via non-equilibrium interacting quan¬ 
tum dot, would require the computation of fluctuation- 
determinant in the vicinity of the instanton solution ob¬ 
tained. Such an analysis is beyond the scope of the 
present paper and is relegated to a separate publication. 
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Appendix A: Korshunov Instantons in Matsubara 
representation 

The dissipative action in the Matsubara representation 
can be written as 

n r 13 

S = - d Tl dT 2 d(r 12 )e^ (Tl) -^ (T2) 

4 In 


+ — I dr (j) — in g 


dr </, 


(Al) 


J c JO 


where j3 = 1/T and t stands for imaginary time. The 
integration kernel here reads 


T 2 

«( T ) = -7T 


1 


+ 


1 


2 [sin (ttTt + id) sin 2 (7rTr — *0) 


■ (A2) 


In order to simplify the comparison with the results ob¬ 
tained in the main text we define the variable 

u = e 2mTT , (A3) 

which is used to rewrite the action (Al) in the form 


5 = 


gT 2 


+ 


2 
zttT 


du i dh e i0(«i)-i0(«2) a ( U i / U2 ) 


2 E, 


o 2niui 2?riu2 
d£ 2 
du 


duu 


c Jo 


du 


d(j) 
du ’ 


(A4) 


where the integration is taken over the unit circle O that 
corresponds to \u\ = 1 and the transformed kernel is 
defined as 


y{u) = 


{u+- l) 2 (u_- l) 2 ’ 


(A5) 


with u± = u( 1 ± <5), S < 1. 

Minimising the action (A4) one derives the saddle 
point equation 


8irT d d(j) 
gE c du du 
1 




,i4>(u')—i(f>(u) _ gi(fr(u)—i(j)(u 


Ll')^ 


(u' — U-) 2 (u' — u + ) 2 


= 0 , 


(A6) 


which can be solved following the original work by 
Korshunov. The solution 0 = <j>i(u), which corresponds 
to an instanton with the winding number one, reads 


e ~ Ci 


u — Z\ 

u-( 2 


1 - 


2r _ 

gE c [(u- C 2) 2 (u — C 1) 2 J 7 

(A7) 


where we neglect the terms of the order of ( T/gE c ) 2 . The 
result of Eq. (A7) is parameterised by a constant c\ and 
complex numbers £ 1,2 such that |£i| >1 and (2 < 1- For 
T/gE c < 1 the solution represents a slightly distorted 
Mobius transform. 

One may further require that the absolute value of the 
bosonic held is fixed hence the phase 4> is rea l- This 
requirement corresponds to further constraint (2 = 1/Ci 
and Ci = (/ 2 e ,A with a real A. Under such a constraint 
the result of Eq. (A7) represents a map of the interior 
of the unit circle to itself. Even though the additional 
constraint on the held seems reasonable we have not been 
able to justify it with a formal argument. 

The leading part of Eq. (A7) describes a circle in the 
complex plane that encompasses the origin just once. 
The sub-leading terms deform the circle without affect¬ 
ing its topology. In the other words, the function still 
corresponds to a closed contour that encompasses the 
origin once. The terms of the order T/gE c are usually 
neglected in literature since they do not affect the value 
of the action at the saddle point solution. 

The action evaluated on the instanton trajectory of 
Eq. (A7) is given by 


Sle^ 1 } 


g 7 t 2 t Ci + C2 

2 4E C Ci — C 2 ’ 


(A8) 


thus reproducing the results of Eqs. (21,22) of the main 
text for the real-time action A = iS. 

Similarly one can construct the instanton trajectory 
(j) = ) with the winding number W larger than one, 


w _ ri 

e i(j>w(u) _ Cw n (A9) 

1 = 1 u ^2 

which is parameterised by an arbitrary number cw and 
2W complex coordinates such that |Cj| > 1, IC 2 I < 1- 
The result of Eq. (A9) corresponds to the action 


S{e i<t,w } 


gW ™ Cj + d 
2 Of 


(A10) 


The leading term in Eq. (A10) is again manifestly equiv¬ 
alent to the result of Eq. (24) of the main text. 


Appendix B: integral transform 

In this appendix we shall summarise some of the in¬ 
tegrals which were used in the main text. An integral 
transform with Cauchy type kernel defined on the con¬ 
tour C leads to the following results 


f dw 1 

< J 

[0, 

n > 0 

(Bla) 

Jc 27r * 

w\ — w — id | 


n < 0 

I 

dw 1 w™ j 

\w n , 

n > 0 

(Bib) 

Jc 

2tt i Wi — w ] 

10, 

n < 0 























7 


For z = z\(w) given by the instanton trajectory of 
Eq. (20) one further obtains 


dw' z 1 (w') 


= 0, 


c 2iri ( w' — w) 2 3 4 5 6 
dw’ z(w') _ 1-& 

c 2ni ( w’ — w) 2 (1 — £w) 2 ’ 

dw' z(w') 


= 0, 


c 2ni ( w ' — w + iS) 2 
dw' z~ 1 (w') 
c 2ni ( w’ — w + iS) 2 (1 — £w) 2 ’ 




(B2a) 

(B2b) 

(B2c) 

(B2d) 


hence one can easily check that the solution ( 20 ) satisfies 
the saddle point equation of Eq. (19). 


Appendix C: saddle-point equation out of 
equilibrium 


The saddle-point equations (19) for the analytic func¬ 
tion z(w), which reduces to z±(w) on the upper (lower) 
side of the brunch cut w £ (—1, 1 ), can be written as 


z 2 (w)I 1 [z](w) = I 2 [z](w), 


(Cl) 


where we defined the integrals 


I\[z]{w) = j) dw' 


1 


+ 


1 


(w — w’) 2 (w — w’ + iSy 


x L 


1 


I 2 [z\{w) = dw' 


1 — w 1 + w'' 

1 4 - w 1 — w’ J z(w') ’ 

1 1 


(C2a) 


+ 


x L 


(w — w') 2 (w — w’ + iS) 2 / 

1-wl + w'\ , 

----- z(w ), (C2b) 

l + wl-w'J K v ' 


where the pole at w' = w lays inside the contour C while 
the pole at w' = w — iS is outside C and the limit <5 —>• 0 
is assumed. The function L{x) is defined in Eq. (29). 
With the help of the Ansatz of Eq. (31) supplemented 
with the assumption on analytic properties of the func¬ 
tion F{w) (such that all zeroes and singularities of F(w) 


are confined to the interior of the contour C) one finds 
' 1 — w 1 + w' \ 1 


h =- 
I 2 = — 2-7T* 


d 
dw 

d 

dw’ 


L 


1 + w 1 — w' J z(w') \ w , 


, (C3a) 


1 — w 1 + w 1 
1 + w 1 — w’ 


7 ) z ( w ') 


J W —YW 


2 il{£ ~ r])F(r)) L f 1 - w 1 + r/ 


(w — rf) 2 


1 + w 1 — 77 


(C3b) 


where we introduced the parameter r] = l/£ for con¬ 
venience. With the help Eqs. (C3) we readily rewrite 
Eq. (Cl) in the trivial form 


dz(w) 

dw 


= r){v ~ 0 F (v)L 


1 — W 1 + Tj 
1 + W 1 — 7] J 


(C4) 


that is solved by Eq. (33). The following indefinite table 
integral has been used for the integration of Eq. (C4) 


dw 


1 — w 1 + 77 
(w — rj) 2 \l + w 1 — r] 


1 — w 1 + 77 
rf — l”" \1 + w 1 — rj 
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const, (C5) 


where we defined 
B v (x) = x l 


1 — x 


— iF 2 (\,iv\ 1 + iv\ x) ) . (C 6 ) 


The result of integration in Eq. (C4) is, then, written as 
v(t~v)F(v) fgl+dR . n9L9R w — r] 


z{w) = C + 
x B v 


w — r/ 
1 — w 1 + rj 


+ 2 - 


g 2 rf — 1 


1 + w 1 — g 


B- v 


1 — w 1 + r) 

1 + w 1 — g J 


■ (C7) 


To obtain the result in the form of Eq. (33) we have 
redefined the integration constant C as 


C = K -77 + 2 


,9L9RV{n~0 


KR , 


g 2 77 - 1 

where K = Fig 7 ) and R is a new integration constant. 


(C 8 ) 
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